{
"cells": [
{
"cell_type": "markdown",
"id": "d26e3b44",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"# Decomposing Qiskit circuits"
]
},
{
"cell_type": "markdown",
"id": "f5994d25",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"In this notebook, we show how Qiskit circuits can be converted into Perceval circuits. To do so, we take the example of a simple gate-based circuit (a circuit producing GHZ states) and we show the translation to a linear optical circuit. We also show the equivalence between the two circuits."
]
},
{
"cell_type": "markdown",
"id": "3ed1c79e",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"As usual we start by imported the useful libraries. Note that this notebook requires the installation of Qiskit (which can be easiliy done with `pip install qiskit`)."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "9dea723b",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"import perceval as pcvl\n",
"from perceval.components import catalog\n",
"from perceval.converters import QiskitConverter\n",
"\n",
"from qiskit import QuantumCircuit\n",
"from qiskit.quantum_info import Statevector"
]
},
{
"cell_type": "markdown",
"id": "61e8fe05",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"## GHZ State generation in Qiskit"
]
},
{
"cell_type": "markdown",
"id": "57423a74",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"We first define the circuit generating GHZ states of 3 qubits with Qiskit. To do so, we first act with a Hadamard gate on qubit 0 to put in superposition of state $|0\\rangle$ and $|1\\rangle$. Then we perform two CNOT gates using qubit 0 as control and qubits 1 and 2 as targets."
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "aef291d3",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/plain": " ┌───┐ \nq_0: ┤ H ├──■────■──\n └───┘┌─┴─┐ │ \nq_1: ─────┤ X ├──┼──\n └───┘┌─┴─┐\nq_2: ──────────┤ X ├\n └───┘",
"text/html": "
┌───┐ \nq_0: ┤ H ├──■────■──\n └───┘┌─┴─┐ │ \nq_1: ─────┤ X ├──┼──\n └───┘┌─┴─┐\nq_2: ──────────┤ X ├\n └───┘
"
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Create a Quantum Circuit acting on the q register\n",
"circuit = QuantumCircuit(3)\n",
"\n",
"# Add a H gate on qubit 0\n",
"circuit.h(0)\n",
"\n",
"# Add CX (CNOT) gates on control qubit 0 and target qubits 1 and 2\n",
"circuit.cx(0, 1)\n",
"circuit.cx(0, 2)\n",
"\n",
"# Draw the circuit\n",
"circuit.draw()"
]
},
{
"cell_type": "markdown",
"id": "cca607ca",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"We display the final state when starting from the input state $|000\\rangle$."
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "7edaa2b3",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/plain": "",
"text/latex": "$$\\frac{\\sqrt{2}}{2} |000\\rangle+\\frac{\\sqrt{2}}{2} |111\\rangle$$"
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Set the initial state of the simulator to the ground state using from_int\n",
"state = Statevector.from_int(0, 2**3)\n",
"\n",
"# Evolve the state by the quantum circuit\n",
"state = state.evolve(circuit)\n",
"\n",
"#draw using latex\n",
"state.draw('latex')"
]
},
{
"cell_type": "markdown",
"id": "9a4a941f",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"## Conversion to Perceval"
]
},
{
"cell_type": "markdown",
"id": "9d6365aa",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"With the use of `QiskitConverter`, we can transform the Qiskit circuit into a Perceval circuit. It uses 2 modes per qubit and additional modes for ancillary photons to perform deterministically two-qubit gates. Below the first six modes correspond to the three logical qubits (see the 'Spatial Modes encoding' paragraph in the 'Basics' section of the documentation) of the gate-based circuit above.\n",
"\n",
"The other modes are used to successfully implement two-qubit gates via heralding or post-selection. Heralding employs [4 ancillary modes](https://doi.org/10.1073/pnas.1018839108) while post-selection employs [2 ancillary modes](https://journals.aps.org/pra/abstract/10.1103/PhysRevA.65.062324). With the option `use_postselection=True` in the method `.convert` on a `QiskitConverter` object, every CNOT but the last is implemented with a heralding scheme. Here it means that it would add $4+2$ ancillary modes. The option `use_postselection=False` only implements heralded CNOTs. Here it would mean $4+4$ ancillary modes. Note: the `use_postselection` option is `True` by default."
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "dccbd8c9",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/plain": "",
"image/svg+xml": "\n"
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"qiskit_converter = QiskitConverter(catalog, backend_name=\"Naive\")\n",
"quantum_processor = qiskit_converter.convert(circuit, use_postselection=True)\n",
"pcvl.pdisplay(quantum_processor, recursive=True)"
]
},
{
"cell_type": "markdown",
"id": "88ff50bd",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"With this converted circuit, we can now check that the resulting state is the same as before the conversion. By default, the input is the logical state $|000\\rangle_L$. Note that where Qiskit displays state in the order $|q_2q_1q_0\\rangle_L$, Perceval uses the reverse order $|q_0q_1q_2\\rangle_L$, but still shown as Fock states. Here, it doesn't change anything since we end with only $|000\\rangle_L$ and $|111\\rangle_L$ states."
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "8ddb1457",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/plain": "",
"text/html": "